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I . INTRODUCTION 


A large number of linear methods exist for estimating 
the lift and other force and moment components of arbitrary wing 
planforms at subsonic speeds. The vortex lattice technique of 
references 1 and 2 and the distributed singularity methods of 
references 3 and 4 have demonstrated both versatility and ac- 
curacy and have contributed significantly to the design process. 
However, most of the analytical work done has been limited to the 
linear or near-linear lift coefficient range and have pertained 
mainly to symmetric flight. Also, the above-named references are 
fairly sophisticated and require a significant amount of calcula- 
tion time on a large computer. 

The present study was aimed more specifically to the require- 
ments of General Aviation aircraft, where the sweepback angles 
are smaller, the airfoil sections conform more towards the NACA 
series of airfoils, the flow is ingompressible, and less computer 
availability is assumed. It was considered worthwhile also to 
provide a capability in the subject computer program for estima- 
ting the aerodynamic force and moment coefficients and stability 
derivatives in asymmetric (yawed) flight and in steady maneuver- 
ing flight as well as in symmetric equilibrium flight. An 
attempt was made also to extend the angle of attack range into 
the nonlinear region near stall. 

The finite-step method described in references 5 and 6 is 
adopted in the subject study. The wing span is divided into a 
number of segments and a single horseshoe vortex is used to 
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simulate each of these segments, in this fashion a fewer number 
of simultaneous equations are required to determine the 
spanwise distribution of circulation than in those methods, . 
References 1-4, where vortices are distributed in the chordwise 
direction- also. It was assumed that sideslipping motion could 
be treated by using horseshoe vortices with 1 cranked 1 trailing 
elements. The trailing elements were assumed to follow the 
chord lines from the line of quarter chords until the trailing 
edge of the wing was reached. Following the trailing edge the 
vortex lines were rotated so that they lie in planes parallel 
to the (yawed) free stream velocity vector. This assumption 
was based on the Weissinger method as described in reference 7. 

It was further assumed that the position and shape of the cranked 
horseshoe vortices were unaffected by the steady angular velocity 
of the wing in maneuvering flight. 

The nonlinear lift estimation in the stall region was to be 
accomplished by replacing the estimated linear circulation 
strength at each spanwise station by the strength corresponding 
to an airfoil section, as determined experimentally, at the 
effective angle of attack calculated at the quarter chord. This 
scheme is essentially that used by Sivells and Neely, 
reference 8. However, the finite-step method was used instead 
of lifting -line theory as in reference 8. 

Several iterations were made to adjust the span loading 
after substitution of the (nonlinear) section data. The process 
was found to be convergent at the .lower angles of attack, but 
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divergent as soon as the effective angle of attack at any span- 
wise station exceeded the stall angle of the airfoil at that 
station. Recent papers dealing with this problem, references 9 
and 10, have shown similar difficulties and reference 9 indicates 
a method for obviating them. 

A computer program has been prepared, based upon the 
above analysis, which has the capability of estimating the inviscid 
incompressible, linear five component aerodynamic coefficients 
of arbitrary wing planforms at a general angle of attack and 
sideslip angle, while undergoing steady rotary motion (P, Q, R) . 

The wing geometry was limited to a single trapezoidal shaped 
panel per side, and includes sweepback, twist, differential 
twist (left- to right-hand wing panels) and dihedral. The wing 
can have arbitrary root and tip airfoil sections. However, 
high lift devices and/or control surface are precluded, as are 
fuselage and nacelle strut interference, and propulsive slip- 
stream effects. 

Numerous computer runs were made on a series of wings of 
different geometry and the stability derivatives calculated were 
compared with the available methods for estimating these deriv- 
atives in the linear angle- of attack range. Very good agreement 
is shown in all but a few derivatives and these are believed 
to be of secondary importance. 

The computer program should prove most useful for straight 
and relatively unswept wings (A<35°) of moderate to high aspect 
ratio (AR>2) . Low aspect ratio wings and highly swept wings 
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will introduce separation at the tips and along the leading 
edges at moderate and high angles of attach, neither of which 
effect is taken into account in the present method. 
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II. SYMBOLS 


A / AR 
A 
b 
c 

C T 

C R 

C MI 

C PL 


m 

i 

C A 

C D 

S, 

c m 

St 

Sr 

c i 

C n 


c/4 


Aspect Ratio, b 2 /S ^ 

ref 

Axial Force 

wing span measured in x,y plane 
local chord line 
tip chord 
root chord 

wing chord coincident with left-hand side of wing strip 
wing chord coincident with right-hand side of wing strip 

Jb/2 

mean aerodynamic chord, 5=2^ o 2 dy/S ref , (reference chord) 

section drag coefficients d */ q c 

^oo 

section lift coefficient (airfoil lift coefficient): 

L */q c 
oo 

section moment coefficient measured about quarter 
chord pc Lnt, (m'/q c 2 ) 

wing axial coefficient: (A /ct s ) 

co ref' 

wing drag coefficient: (D/cr s ) 

' co ref ‘ 

wing lift coefficient: (L/q S ) 

oo ref ; 

wing moment coefficient: (m/cr s c) 

u/4 co ;3 ref CJ 

wing normal force coefficient: (N/a s 

• ' ^co ref 

side force coefficient: (Y/q S f b) 


) 


rolling moment coefficient: (z/q S 


'co ref 


Jb) 


yawing moment coefficient: (n/cr s b} 

v /y oo a ref j 

Sf 9 S i point°Sn a c/4 SPlaCement vector ' or assignation oi 

S^SivrSinn^to? aSmea parallal fc ° stream 


O^IGJNAL sage IE 

OF POOR QUALITY 






SP 

p 


f/g*h 


i 

A £ A 

j 

k l 

K 

^ref 

L 


m 

M 


M 

n 

N 


influence coefficients used for calculating the 
induced velocity 
elementary force vector 
force vector 

coordinates of a point measured with respect to the 
origin of a horseshoe vortex, see figure 6 
normal distance of a point measured to the line of 
action of a vortex segment 

an integer used to count horseshoe vortices, from 
left to right 

unit vector triad associated with x,y,z axes 
an integer used to count wing strips 
factor of proportionality for reducing angle of 
attack: O^k^^l, see Appendix A 

an integer used to count the segments of a horseshoe 
vortex 

vortex segment length of trailing elements 
general reference length 

total lift force, measured perpendicular to free 
stream relative wind vector 

total pitching moment measured in plane of symmetry 

moment vector measured relative to position of the 

center of gravity 

Mach number 

total yawing moment 

number of spanwise strip (equals horseshoe vortices) 
considered. Value includes strips on both left- and 
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A 

n 3/4 


P,Q,R 


AJ Arf AJ 

P,Q,R 


q 


CO 





A 

r l/4 


A 

r 3/4 



r i,r 3 



right-hand sides of wing 
normal force component 

unit vector perpendicular to local wing chord plane 
with components n^ , n^, n^ on left, and n^ g , n 12 

on right-hand sides of wing 

unit vector normal to mean camber surface of wing at 
the 3/4 chord line with components n 1# n 2 , n 3 on 
left, and n ? , n Q , n g on right-hand side of wing 
rolling rate, pitching rate and yawing rate, 
respectively, radians per second 

<NJ \ .... 

nondimens ional rolling rate (P- ^ — ) , pitching 

_ — CO “p 

rate (Q= ) , and yawing rate {R = ~ — ), respectively 

co 03 

free stream dynamic pressure 

velocity vector, magnitude, and unit vector induced 
at a point P owing to a single vortex segment 
unit vector coincident with line of quarter chords 
with components r 4 , r 5 » r g on left, and r 1Q , r.^, 
r 12 on right-hand side of wing 

unit vector coincident with line of three-quarter 

chords with components r^, r^, r 3 on left, and r^, 

r_, r n on right-hand sides of wing 
o y 

displacement vector from' airplane center of gravity 
to the control point "c^" 

displacement vectors from the initial and final ends, 
respectively, of a vortex segment to the point P 
unit vector coincident with a vortex segment, positive 
sign is chosen in agreement with right-hand rule for 
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s 


ref 


A 

H/4 



u, v,w 


v 


U,V,W 


V 




x,y,z 


X,Y,3 


direction of circulation, I\ 

semi-width of wing strip, also semi-span of horseshoe 
vortex, measured in plane of vortex 

reference wing area, the wing planform area projected 
on to x,y plane 

unit vector tangent to the mean camber surface of 
wing measured in a chord plane. The components of 

^■ 2/4 are t 2' fc 3 ° n le:Et ' and t 7' fc a' on 
/ 

hand side of wing 

unit vector coincident with local chord lines with 
components t^_, tg, tg on left, and t^g, t.^, t ^ 
on right-hand side of wing 

induced flow velocity components along the x,y,z axes 
induced velocity at point P due to the entire system of 
horseshoe vortices 

total flow velocity components along the x,y,z axes 
total relative velocity vector existing at a point P 
owing to the sum of the free stream, rigid body 
rotation, and induced velocities 

magnitude and velocity vector of the free stream 
orthogonal axes to which wing geometry is referred. 

The x axis is coincident with the root chord line, 
the y axis is normal to the plane of symmetry of tie 
wing, positive to the right, and z axis is defined 
according to the right-hand rule 
body force components 
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GREEK SYMBOLS 


a 

P 

r 



fi 


G 

0 

X 

A 


P 

cp 


. * 


(JU 


angle of attack 
angle of sideslip 

dihedral angle, positive for wing tip raised above the 
xy reference plane. 

vortex strength of the x horseshoe vortex. 

= r./v m 

X 03 

angle measured between mean camber line and chord line, 
positive for positive camber (concave side on bottom of 
wing) . 

local twist angle 

acute angle between and r 2 or obtuse angle between 
R 3 and r 2 * 

taper ratio of wing ( 1 = c^/ c R ) 

angle of sweepback of a constant percent chord line. 

The quarter chord line is designated without a subscript? 
other lines are subscripted according to fraction of 
chord , 

atmospheric density 

angle of inclination when viewed from a line parallel 
to the x axis between bound vortex segment and y axis . 
Positive for clockwise rotation when looking forward. 

angle of yaw of bound vortex segment. Positive for 
counter clockwise rotation about z axis when viewed from 
above. 

angular rate of rotation of wing in space, components are 
F, Q, R. 
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SUBSCRIPTS 


( >B 

body axes 


(, v 

control points located at 3/4 chord line at midspan 
the wing strips. 

^ ^ CG 

center of gravity 

of airplane. 

< ) d . 
J 

< > e 

m 

points coincident 

with origin of horseshoe vortices 

points located at 
edges of the wing 

the 5/8 chord line at the lateral 
strips . 

< >i 

denotes the i horseshoe vortex. 

■ ( >3 

l.'L 

denotes the j spanwise strip. 

< >1 

local value. 


( >1 

left hand side. 


< )» 

an integer used to 

count chordwise vortex segments. 

< >o 

root chord value 


* ^R 

right hand side. 


^ ^ RB 

rigid body value. 


^ ^S 

stability axes. 


( ) T 

tip chord value. 



wind axes. 


SUPERSCRIPTS 


o 

a vector quantity. 

(except for c) . 

( A ) 

a unit vector. 


()' 

two dimensional or 

section value. 

•ABBREVIATIONS 


L.E . 

leading, edge. 


T.E . 

trailing edge. 
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III. ANALYSIS 


The various calculations which were required for the 
preparation of the computer program are presented in this section. 
These calculations have been broken down into three main cate- 
gories, namely: (a) the determination of the wing geometrical 

characteristics as required for determining the "no-flow" 
boundary conditions, and calculating the effective angle of attack 
as required in the nonlinear analysis, (b) the finite- step linear 
analysis, and (c) the modified finite-step method which incorpor- 
ates the nonlinear airfoil section data. 

A. Wing Geometry 

Figure 1 presents the planform geometry of the wing and 
establishes the body axis system used. The x,y plane lies in the 
plane of the page and hence only the projection of the free stream 
velocity is shown. The wing shown in figure 1 is planar in that 
neither twist nor dihedral are indicated. 

, Figure 2 presents a perspective drawing of the left hand 

panel of the wing which has been twisted, (zero dihedral case) . 

It has been assumed in this figure: (a) that the root chord lies 

along the x axis, and (b) that twist is introduced by rotating the 
tip chord in a plane parallel to the plane of symmetry about the 
quarter chord point. The tip chord twist angle is defined as e T , 
and is subscripted with an L or R which designates the left- or 
right-hand panels-, respectively. The twist angles of the left or 
right hand panels can have different values in order to study the 
effect of twist asymmetry. Notice that the twist angle is defined 
as being positive for the case of washin. This angle will normally 
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be negative as most wings are washed out. The leading and trail- 
ing edges of the wings are assumed to be straight lines, which 
causes the spanwise variation of the twist angle to be nonlinear 
if the wing is tapered. 

The dihedral angle of the wing is produced by rotating 
uhe quarter chord line (and hence the entire left or right hand 
wing panel) about the x-axis through the angle r. The twist and 
dihedral angles are assumed to be relatively small so that the 
order in which the rotations take place is inconsequential. In 
the analysis which follows the twist angle .rotation is assumed 
to occur prior to the dihedral angle rotation. 

The purpose of this section is to present the results of 
an analysis for (a) calculating the spanwise variation of the 
twist angle, (b) determining the unit tangential and normal vectors 
required to satisfy the "no-flow'* boundary conditions at the three- 
quarter chord control points, and (c) determining the local angles 
of attack and sideslip at the quarter chord line. The left hand 
side of the wing is treated first, and then the right. 

Left hand side of wing , y<0 

The local twist angle of the left hand wing panel, e^, at 
the y station, (y£0) , is calculated for a wing without dihedral 
using the relationship: 



Z LE Z TE 
a LE~ X TE 


(1) 


Upon determining expressions for the x,z positions of the leading 
and trailing edges as a function of y, and using the definition 
for taper ratio, namely that X. = c^/c^, equation 1 can be expressed 
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\ (sin e ) 


■ l L b/2 

tan c T = - vv - 7 s — v — ~ — tLLZ- 
L 1+ (1-A cose_ ) 


? (y^o) 


L ' b/2 


This equation could serve for the right hand panel also, (with 

certain sign changes) , except that in the general case e can 

T L 

differ from e . 

The "no flow" boundary condition requires that the (total) 

flow velocity vector be tangent to the mean camber surface of 

the wing at the three-quarter chord location of the control 

points. This boundary condition is prescribed most easily by 

setting the scalar product of the velocity and local unit normal 

vector, defined as ^ 2 / 4 ’ equal to zero. In order to calculate 

the unit normal vector, the 3/4 chord line unit vector, r_ , 

3/4 

and the local unit tangent vector, 1 3/ , 4 , must be determined. 

Figure 3 presents a schematic . drawing of the vector triad, namely 

r 3/4 ' t 3/4 and n 3/4 which is used in the analysis given below. 

L L L 

Figure 4 shows a cambered airfoil section at a typical 
spanwise station prior to the dihedral angle rotation. The inclina- 
tion of the mean camber line relative to the chord line is defined 
as ^3/4 anc * ^is an Sfle is assumed to have a linear spanwise varia- 
tion as given by: 

, , . i„ l ORIGINAL PAGE lt> 

S 3/4 = V4 0 "[ 6 3/4 0 _6 3/4 t ]e/2 OF POOR QUALITY (3) 

where & 3 /4 is tlie inclination measured at the 3/4 chord station 
of the root airfoil section and 5 2/4 ^- s the value measured at the 

tip airfoil section. Equation 3 is seen to apply to both the left 
and right hand wing panels. If the zero lift angle of the root 
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and tip chord airfoil sections are Known, these are substituted 


for 5 3 y 4 and « 3/ , 4 , respectively. 

The unit tangential vector, is rotated about the x 

axis through the dihedral angle r to bring it to its final ori- 
entation, namely . This latter vector is defined as: 


A 

t 


3/4 = *1? + t 2^ + fc 3* 
L 


which may be expressed more compactly as: 


where 


fc 3/4L “ H' t 2' t 3 


t x = —cos (e L + 5 3 / 4 ) 
t 2 = -sin (E L +5 3 ^/ 4 )sin r 
t 3 = sin (^+63^)003 r 


(4) 


The displacement vector R 


3/4 t 


is defined as that vector 


whose tail lies on the 3/4 chord point of the root chord and whose 
head lies on the 3/4 chord point of the left tip chord. It may 
be considered as that vector which passes through the locus of 
control points and its magnitude is given by: 

• ^ =[(? - I tan * - T- COSe T ) + (!) + (? sine T - I tan ] <Sa) 


The vector r 


3/ 4 t 


“L “ L 

is the unit vector associated with the displace- 


ment vector ’ 8 - 2 / 4 : an< ^ as: 

L p 

- ' . 3/4 l _ 


•3/4. 


L 


| 3/4_ ( 


= r r r 2 , r 3 


(5b) 
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where 


r i =( T 1 ~ I tan A " 5 s cos 

b/2 

T* *~~* r i r 


(5c) 

(53) 


r 3 ” ( 2^ sin e T - tan r)/3j 

A 


(5e) 


The -unit normal vector, n 3 y 4 , is calculated as the vector 

A A 

product of r„ and t_ and is positive in the upwards direc- 
tion, (see Figure 3) : 


A A 

— r 


A / A A 

x t- /A ./ .r n x t 


n 3/4_ X 3/4 t x c 3/ 4^' / j~3/4 r 

ij J-l JLI i-1 


3/4 t 


(6a). 


and is defined as: 


a _ 

n 3/4 T " n l' n 2 x n 3 
Xj 


(6b) 


where 


n l ■ (r 2 t 3 " r 3 t 2 )/d 3/4. 


L 


(6c) 


n 2 ^3^1 " r l t 3 ^ / d 3/4 T 


(63) 


n 3 ^ r l fc 2 ~ r 2 t l^ //d 3/4 


(6e) 


and 

d 3/4 = ^ r 2 t 3~ r 3 t 2^ + ^ r 3 t l“ r l t 3^ 

L 

+ ^ r l t 2“ r 2’ t l^ ^ 2 “ 

In order to calculate the "effective" angle of attach, of 
this airfoil at the quarter chord point as required by the non- 
linear analysis, the vector triad associated with the quarter 
chord point r^y^ , , and n^y^ . is calculated. This triad 


(6f ) 
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has the same general orientation as the three-quarter chord triad 
which is shown in Figure 3. Note that the unit tangent vector, 
t l/4 L ' is now tan ? ent to the chord line, whereas ^ 3/4 was tan- 
gent to the mean ' camber line. The quarter chord triad is defined 
as; 


r l/4 L s r 4' r 5 ' r 6 

(7) 

^/\ S ' ±5 ' 

(8) 

n l/4 L = n 4' n 5' n 6 

(9) 

r 4 “ -tanA/Cl+tan 2 A+tan 2 r] ?5 

(10a) 

r 5 = -l/Cl+tan 2 A+tan 2 r]^ 

(10b) 

r 6 = ~tanr/Cl+tan 2 A+tan 2 rT* 

(10c) 

t 4 = -cos Gl 

(11) 

t 5 = -sin e L sin r 

(11b) 

t 6 ~ sin ^ cos r 

(11c) 

n 4 ", <^'-6 - Vs^lAj. 

(12a) 

K 5 “ <*6*4 ~ r 4'T,>/'V4 L 

(12b) 

n 6 = (% - r S t 4 )/^ 1/4 

L 

(12c) 


and 

VV C(r 5 t 6 " r 6 t 5 )2+(r 6V r 4 t 6 )2+ < r 4 t 5- r 5 t 4 ) 2 ] % 
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Ricrht hand side of wing, yaO 

The subscript ( ) will be used below to designate the 
right hand panel as indicated earlier. The equations to be 
given reflect the mirror image of the left and right hand panels 
except for the fact that the wing tip twist angles can have dif- 
ferent values ■ The development of the equations for the right 
wing panel parallels the analysis given above. 

The local twist angle for the right hand wing panel is 

given by: 

tan 

which, when compared with equation 2, shows sign and angle changes. 

The unit vectors associated with the vector triads located 
at the three- and one-quarter chord stations follow directly: 

t 3/4 _ s * 7 * t 8' fc 9 
R 


Msin e^) 

g R = l-[l-A.cos e ^ /2 


7 (y^o) 


(13) 


where 


t ? = -cos ( e R + S 3/4 ) 


tg = sin (e R +fi 3 / 4 ) sin r 


t g = sin U r +6 3 / 4 ) cos T 

The magnitude of the displacement vector Rg^^ i s 9'^ ven ^y 

R 

d* = f(|t B nA+. ^ cos ^ - 3&) 2 + (§)+.(§ tan T 
*- R 


- sin s ) 
^ X R 






(14a) 

(14b) 

(14c) 


(IS) 
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The unit three-quarter chord vector r 3/4 of the. right hand panel 
is defined as: 


r 3/4 R “ r 7' -V r 9 

Where 

„ _ ,°R b g t . 

r 7 ~ ^ - 2 tan A ~ ~ COS e T )/d R 

R 

r 8 - ! 

c ijt Id 

r g = ( 2 “ sin R T ~ 2 tan r ^/ d R 

R 

From Figure 5 it is seen that: 

= V 4 j, X S 3 / 4 r 

— * p-j 8 9 

where 

n 7 = (t 8 r 9 “ t 9 r 8 )//d 3/4 

R 

n 8 “ ^ 9 X 7 ~ t 7 r 9^ //d 3/4 

R 

n g = (t 7 r Q - t 8 r ? )/d 3/ 

R 

and 

d 3/4 R ~ ^ t 8 r 9“' t 9 r 8^ + ^ t 9 r 7“ t 7 r 9) 2 
+ (t 7 r 8- t 8 r 7 )2]% 

The vector triad associated with the quarter-chord point: 
. ^*1/ 4 r s t 10 / t ll / t 12 

where 

t 10 = -cos 


(16a) 

( 16 b) 

(16c) 

(17a) 

(17b) 

(17c) 

(173) 

(17e) 


(18a) 


18 



where 


and 


where 


t, . = sin sin 

JLX R 

r 

(18h) 

t X2 = sin e R cos 

r 

(18c) 

*l/4 R 55 r 10' r ll 

' r 12 


r 10 = r 4 


(19a) 

r ll = " r 5 


(19b) 

r 12 = r 6 

P 

(19c) 

*±/ 4 ^~ n 10' n ll' 

R 

n 12 = t l/4 R X X 1/4 r 


n 10 “ (t ll r l2 

t 12 r ll ) / ’ d l/4 R 

(20a) 

n ll “ (t 12 r 10 

t 10 r !2 )/d l/4 R 

(20b) 

n 12 = (t 10 r ll 


(20c) 


d l/4 “ £ (t ll r 12~ t 12 r ll ) ' 1 ' (t 12 r 10" t 10 r 12 ) 
R 


2 -M 




(20d) 


Sufficient geometric information has' heen determined so that 
the finite-step method and its subsequent non-linear modification 
can he applied. 
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Be Finite-Step Method, Linear Case 

The spanwise variation of circulation along, a wing is 
estimated in the Finite-Step Method using a set of single horseshoe 
vortices in the chordwise direction and H such sets of vortices dis- 
tributed in the spanwise direction, a reasonable number for H 
is 40, which corresponds to 20 vortices per side of the wing. 

The finite-step calculation used here is similar to that pre- 
sented by Campbell (reference 5) and Blackwell (reference 6) . 

The horseshoe vortex elements considered herein have bound 
elements which are both swept and inclined to the x,y plane as 
those of Margason and Lamar (reference 1) , but lie in the pi 
of the. wing rather than being centered in the x,y plane of the 
wing as assumed in reference (1) . The trailing legs of the 
horseshoe vortices follow the chord lines initially and then 
have a crank at the trailing edge of the wing, see Figure 6 . 

Those segments which are located downstream of the trailing edge 
lie more nearly in the free stream direction. In the case of 
sideslip these segments have a deflection angle p with respect 
to the plane of symmetry, so they are streamwise with respect 
to the sideslip variation. In the case of angle of attack the 
trailing segments have a deflection which is proportional to the 
angle of attack in order to simulate the Kutta-Joukowski condi- 
tion at the trailing edge. (It has been suggested that the 
constant of proportionality be one-half in reference 2.) However, 
the constant of proportionality, ^ , can be varied at the option 
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of the user.* 

Figure 6 presents the plan view of the wing and illustrates 

the spanwise distribution of the horseshoe vortices and shows a 

typical wing segment. The semi-width of the horseshoe vortex 

is defined as s, where s is calculated: 

s = — ^ (21) 

S 2Ncosr v 

Referring to this figure, it can be seen that the point x.^, 
a. is the origin of the i th vortex. This vortex has a spanwise 
or bound segment which is coincident with the line of quarter 
chords of the wing and is therefore displaced vertically from 
the horizontal reference plane, the x,y plane. The axes of the 
i^k vortex, namely, f,g,h, are parallel to the x,y,z axes but 
are displaced by the distances x^, The chordwise seg- 

ments are assumed to be parallel to the x axis , and hence the 
twist angle is not simulated exactly. The maximum twist angle 
of the typical wing is less than 5° everywhere, so that the 

physical displacement of the horseshoe vortex and wing segment 

. . , , .th 

should not be significant. The boundary condition on the j 

segment at the control point x c ., y c . , z has been moved to the 

^ 3 • A 

plane of the vortex. (The inclination of the unit vector n^y^, 
which defines the normal direction to the mean camber surface or 
the wing, has been calculated in the preceding section. The 

*The value 1c. equals zero causes the downstream segments to lie 
in the plane of the wing, and the value 3c a equals unity causes 
these segments to be parallel to the free-stream velocity vector. 
Further details are given in Appendix A. 
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components of & 3/4 are determined at the control point which is 
designated with the subscript { ) ) «. 

*he velocity induced at the j th control point due to the 
i th horseshoe vortex is calculated according to an algorithm 
Which is presented in Appendix A. The general horseshoe vortex 
is considered in this Appendix to have five straight-line seg- 
ments: a bound or spanwise segment which was previously described 

as lying along the quarter-chord line of the wing, two chordwise 
segments which lie on, (or near, if the wing has twist), the left 
and right-hand edges of each segment, and two segments which 
extend downstream from the trailing edge to a distance of 
approximately one-thousand root chords, which approximates 
infinity. These five segments are designated K=l-*5 where the 
segments are counted clockwise when viewed from above. They will 
be referred to as the bound, chordwise, or trailing segments in 
the discussions which follow. 


The velocity induced at the point P by the i 
vortex can be expressed functionally in the form 




horseshoe 


_ h 


[, 

L u ■ 


P(f,g,h) 4tr (p) i+P V£^ p j 


j+F. 




w 


r i, (P) 


where the influence coefficients are defined as F u , F v , and F 

These influence coefficients can, in turn, be expressed in 
fxinctional form: 

F=F (f ,g,h, s, cp, a , c MI ,c pL ) 

The analytical egressions of these functions will be exceedingly 
lengthy and while impressive, have no particular virtue in them- 
selves. Therefore the velocity induced by each vortex segment 
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is calculated separately and then summed by the computer for 

For the special case of p=0 and k a =0, the calculations 
agree with results obtained using the influence coefficient 
formulas given in references (1) , (5) and (6) . 

The spanwise variation of circulation is determined by 
calculating the induced velocity at a control point, "c^" , for 
the N system of horseshoe vortices where each of these vortices 
has an assigned but unknown strength It. This velocity is then 
added to the sum of the free stream velocity and the velocity 
produced by the angular rotation vector m at c - The boundary 
condition requires that the total velocity vector is tangent 
to the mean camber surface at the control point. The control 
point is located at the midspan, three-quarter-chord line of each 
wing strip. The simultaneous matching of the flow tangency re- 
quirement at the N control points leads to a set of N linear alge- 
braic equations which are solved for the vortex strength values. 

The Kutta-Joukowski theorem is applied to calculate the incremental 
forces which act on the bound and chordwise segments of the 
horseshoe vortices. Incremental moments due to these elementary 
forces are determined also using the center of gravity as refer- 
ence center. The incremental forces and moments are integrated 
over the entire wing in order to determine the force and moment 
coefficients desired for both body and stability axes. 

1* Calculation of the Induced Velocity, v c _. 

til 

The location of the origin of the i vortex, counting from 
i=l on the left wing tip to i=N on the right wing tip, is given 
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by the following equations: 

Y i = ~ 2 + ( 2 i -l)s cosr 7 i=l to N 

- jyjJtanA 
z i = ~ \y± |tanr 

0 ?he location of the j th control point, » Cj », 


following equations: 


Y Cj 2 4 - ^ 2 3“1) s cosT ? j=l to N 

x =5 = " |y °j l tanA - I CrI 1 - £ 7 )^ ll 


z Cj [y C j [tanr 


The displacement distances are calculated as: 


" " X =j “ X i 
^ = y Cj - y,- 


h ~ z 


z , 


The induced velocity at the control point, "< 

system of i«l to N vortices distributed even] 
span is calculated as: 


where 


% U c • ' V C ' W o 
3 C D G j c j 


N r 

V i p u i:j 

J I - — 1 X J 


i=l 

N 


v. 


= V 
'■j 4 ’Sj 


(22) 

is given by the 


(23) 


(24) 

11 , owing to the 
over the wing 

(25) 
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w. 


K 

b r i 

! j~ L 4 TT 

X=1 


F 


w 


ID 


(26) 


The influence coefficients are given by equations (A24) , 
(A25) , and (A26) , which are presented in Appendix A. 


2. Determination of Spanwise Variation of Circulation 

The unit normal vectors located at the j control point 
have been calculated in the geometry section as: 


S 3/4 l = 

(n 1# n 2 , n a ) 

y Cj < 0 


“3/V 

(n>7 , ^ 9 ) 

y 0j > 0 

(27) 


The velocity at the control point will be tangent to the wing 
surface when the scalar products of the unit normal vectors and 


the total velocity vectors are zero: 

V £ 3 /4 - ° 


(28) 


where the total velocity vector V„ is the sum of the free 

3 

stream relative velocity vector, (V^ ) , the relative velocity 
produced by the angular rate of rotation (-rnXR } , and the 


induced velocity (v c> }: 


V V ®'“”®CCL + % 


(29) 


The boundary condition, equation (28), becomes, after substi- 
tution of equation (29) and rearranging: 

(30) 


ft 3/4 = ( - V co +l,Ki W. > ‘ "3/4 
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The left-hand side of equation (30) contains the unknowns 1^, 
i=l to N, while the right-hand side is known. The velocity on 
the right-hand side of equation (30) is recognized as the 
rigid body velocity at the control point and for simplicity 


the terms are combined: 

Y = - V + wxiL- 
V RB„. CO CG c 


(31) 


where 


A . A .A 

Y = — Y (cos ficosa l+smp j+cospsma k) 
oo co 

«) = Pf Qr R 

5 CG - ' X c 1 - X CG )4+(y C- i - y CG , ^ +1Zc i“ ZcG)fe 
CA J J J 


(32a) 

(32b) 

(33) 


Equations (29) to (33) are normalized by division of the magni- 
tude of the free stream velocity vector. The nondimens ional 
rolling, pitching and yawing rates P, Q, and R are introduced 

also where: 


= , * = J*S- . R = 


Kb 


p “ * Q ” 2V 7 R 2V 


(34) 


2V 


oo 

A- 


co oo 

The nondimensional rigid body velocity components become: 


U. 


RB, 


- 


Z c~ Z CG y c^“ y CG 


oo 


=cosecosa+Q( ) ” R ( 

c/2 V2 


V. 


BE. 


'3 _ 


X c,“ X CG „ Z c j “ Z CG 

J v T*.r 


v 


oo 


=sinp + K (' b/ 2 ) T p ( b/2 5 


(35) 


W. 


RB 


y c-T y CG ^ X Cj“ X CG 


Y^ 2 =cospsina + P( b /~^ “ Q( - /2 } 
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for the 


The left-hand side of equation (30) , divided by , 
left-hand wing panel can be written as: 


v 


N 


V 


oo 


3 A 

• »3/y 


4,1 3=1 to 


i=l 


(36) 


The simultaneous equations for I\ , i— 1 to N are obtained by 
substitution of equations (35) and (36) into equations (30): 



for j= | +1 to iJ 


(37a) 


(37b) 


Equations (37) represents N linear, algebraic equations which 

are used to determine the N values of (T.j./V^) ■ simplicity 

r. is defined as the ratio of T. to V , and this variable 
x x oo 

will be used subsequently. 


3. Calculation of Linear Force and Moment Coefficients 

The elementary force produced by a vortex filament of 
length i, immersed in a flow is calculated using the Kutta— 
Joukowsbi Theorem, as applied by Prandtl: 
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(38) 


5F “ pV XT.G 

■where V is the local velocity vector determined at the mid- . 
k . 

point of the vortex segment under consideration. Equation (38) 
will he used to calculate the forces acting on the ohordwise 
and hound segments separately, in that order. 

a) Chordwise segments 

There are N+l chordwise segments and the subscript m will 

he used to number these segments, see Figure 7. The right- 

iih. . . . . 

hand chordwise segment on the j wing strip is corncxdent 

. th 

with the left-hand chordwise segment on the (j +1) strip. 

The force on the combined segment will he calculated using 

the net circulation acting. The first and last segments, on 

the left and right-hand wing tips, respectively, will have a 

circulation which corresponds to the value existing on the 

horseshoe vortices located at the respective wine tips. The 

chordwise vortex segments are counted using the integer 

variable m, where m varies from 1 to N+l, and the counting 

again proceeds from left to right beginning at the left wing 

th 

tip chord. The net vortex strength acting on the m segment 
are calculated according to the equations given in the table: 


Value 

of m r(m) Value of i 

T 

r m=l _r i“l i=m=l 


m, m=2 -* N 

r =r.-r- , 

m x x— 1 

i=ra 

(39) 

N+l 

r m-N+l” r i=N 

i=N 
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The midpoint of the chordwise segments are subscripted ( ) 
and the induced velocity at these points are* calculated using 
the method described in Appendix A. The induced velocity due 
to a (single) segment at a point which lies anywhere on that 
segment, or on the straight line extension of that segment is 
equal to zero. A test is provided by calculating the perpen- 
dicular distance between every line segment and the point in 
question. If the perpendicular distance is zero (the practical 
criterion used is A-ksO .Is) , the induced velocity from that 
segment Is not calculated, which avoids the indeterminate form 
which would result were the computer to evaluate the equations 
for the induced velocity at that point. 

The position of the points e m are calculated: 

y - “ 2 +2 (m-1) s cosr, m=l to (u+l) 
fc m 

’V - 1 c k-ly era l tan Vs (40) 

z e = - |y e |tanr 
E m m 

where 

ta *v 8 =tanA - ■ (4i) 


The induced velocity at the point e m is calculated: 
v N 

e m _ 1 V ? r 4 A a 

4tt L r i [F u 1+F v I]+F w k] 
i-1 


using: 

f ~ x -x. 

3 . 
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(42) 


( 43 ) 




where it are those values obtained by the solution of equation 
(37). The influence coefficients are again those defined by 
equations (A24) { (&2E), and (A2Q. 

The vector r m is determined using the magnitudes T m 

A A 

as r>er equation (39), and the unit vectors r„ =i which correspond 

K=2 

to the left-hand chordwise vortex segments (K=2) as defined in 
Appendix A, except for the right-hand wing tip segment (m=Ntl) 
where the negative of the 'i unit vector is used: 


•^m *m ^ ' 


m = 1 to N 


and 


(44) 


%fl “ r N+l 1 


. m = N + 1 


The lengths of the chordwise segments are calculated 


and 


^m 4 C MI . 


Vf-l 4 (C PL } n 


m = 1 to N, m = j 


4 C T 


(45) 


The total velocity is calculated at the points e 


m 


V = v - v 
e_ e„ RB 
mm e, 


m = 1 to N + 1 


m 


(46), 


where is defined by equation (3l) except that the point of 

e m 

determination is now e rather than c- as in that equation. The 

ra 3 

terms or components of V__ /V are given by equation (35). upon 

e m 

replacing the .values x . y , z with x , y , z 

c 4 c/ c . e J e ' e m 

i i i m m m 


The incremental forces on the chordwise segments are 
calculated by substitution of the results of equations (39). , (44) , 
(45) and (46) into equation (38) , to wit: 
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(■ 


SF\ = 


D ' e_ *~' k m 

K m 


V 3 civ i) ('T ® ! _ ) * m— 1 to N 

v - xu m ; v 4 Ml'm. 
m 


and 


( 4 ) - 

\D 'wil 


v x(r N i) | c T ? m=N+1 
■p 'n+i ®s+i N 


(47a) 


(47b) 


These expressions may be nondimensionalized by division of both 


sides of the equation by | S ref- . Defining the free-stream 
dynamic pressure, q^ = \ pV^ , the nondimensional incremental 
force equation becomes: 


V 


( aF \ _ J 2 (3 \ \ — m x i ; m=l to N ( 48 a) 

Vq^ J " S -\4 C Miy m \ i m/V m 

'Vref m ref m “ 


and 


V 

■f SB \ _ _2 (3 „ Wr \ gg+ - 1 - x i : m=N+l (48b) 

Wrefl., S ref U * M V - 


N+l 


The total force produced by the chordwise segments is: 

N __ 

p ^ f N . V / 5F_ 


^ q ro S ref^ 4 ^ m S ref 


CH 


‘.’ref' , ' q => S ref ‘ 

N+l m=l 


(48c). 


The resolution of the body axes force components along stability 
axes is given following the derivation of the bound segment force 

vectors . 

The incremental moments, referenced to the center of 
gravity, produced by the chordwise segment is: 

«V = 5 CG x (49a ’ 

e m 

The body axes moment coefficients are defined by the equation: 

“ = % S r= { i + S 3 - f 9 " Cref ' ^ + ^ q " b ( l b) 
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The moment is nondimens ionalized by the factor q oo S ref -t ref 
and is calculated using equations (48c) and (49) : 


M 


^ g co S ref^ref> 


CH 


CG 


b N+l 


^ref 


x( 


6F. 


N+l 


^co^ref 


N R CG. 


m 


m=l 


^ref 


x( 


6F 


m 


^•co S ref 


) 


(50) 


b) Bound vortex segments 

There are elemental forces on the N bound vortex segments 
which must be summed in order to find their contributions to 
the overall forces and moments. The increment forces are ob- 
tained using equation (38 ) , where the local velocity is again 
evaluated at the midpoint of the segment, which in this case 
is the origin of the horseshoe vortices. This point is defined 
as point d, and the coordinates of d are given by equation (22). 
For the sake of completeness, the displacements and the induced 
velocity at the point d are written: 


f = x i 


* “ y i 


h “ z clT z i 

3 


where 


(51a) 


ORIGINAL BAGE 19 
OF POOR QUALrHI 


y d = - 2 

3 

x d = - |y d > anA 

3 3 


+ (2 j— 1) s cosr , y.= - \ +(2i~l)s cosr 




= - Iy d |tanr 
j 

j=l to N 


x 2 
x ± = - |y^ | tan A 

Z ± = - [y^ |tanF 

i=l to N 


(51b) 


and the induced velocity is calculated: 
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■is X f iD r «* +p v* +F w fi ] 


where F u , F v , F w are summed over K as shown in equations (A 24) , 

(A25) , and (A26) , except the terms K=3 are omitted, since the 

point dj lies on hound vortex segments a The total velocity is 

calculated at the point d . ; 

3 


V V VRB d. 


(52) 


where V is defined by equation (31) except that the point of 

a j 

determination is now d^ rather than c^ as in that equation. The 
components of V RBd, A ro are given by equation (35) upon re- 
placing the values (x,y,z) c with (x,y,z) rf as given by equation 

3 j 

(51b) . 

The length of the spanwise segments, and their circula- 
tion strengths are given by; 

■l = 2s [l+cos 3 Ptan 3 (53) 


r j r i rH j,K=3 ? 


These results, equations (52), (53), and (54), are substituted 
into (38) .and the increment force on the bound vortex 
segment is calculated; 


{ '7 i) BND = x ij (2s [l+cos 2 rtan 3 A]^) ; j«l to N 
Division by ^ V cQ S ref leads to the nondimens ional equation; 
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4s 


5F. 

(— 3 — ) = . 45 , 

„ oo ref. ref 


v d _ 

[l+cos s rtan 2 A*P(r.' .) r — -l- x r~ 

' J 1 3“i J V * ra j / K=3 
- CO J 


. . j=l to N 

The total force on the hound segments is given by: 

• N 

-) -i / j=i 


(55) 


(■ — — — ) = Y f — 

~ e = B ® L ( qco s ref ' 3 


q s _ 
^•oo ref 


j=l 


(56) 


The incremental moment coefficient, measured about the c.g. 
of the j bound segment is: 

R 




CG . 


6F- 


U - L " "1 

(~ 5 7 ) ■nvm - "7 7C (“ ~ ) 


9oo s re*W Bm > *re£ 


■Wref BIro 


(57) 


The total moment is obtained by summing equation (57) : 

N 


( « ) = y f 

g oo S refW BND A 


6Mj 






co ref vref 


(58) 


Finally we sum the force and moment coefficients on the 
chordwise and bound segments using the results given by (48c) 
and (56) , and (50) and (58) , respectively: 


(. F ) = ( I \ j_ , F 

n / m \, 


g co S ref' T ^BSiD + { q S 


-) 


l co u ref avilJ ^co Q ref CH 


(59a) 


M 




M 


^refW 1- V 5oo S rertre£ )BHr,+ ( «co s refW ,CH 


M 


(59b) 


The total force is defined using body axes components 


- /Jv A A 

F T = Xi + Yj + zk 


where X is the negative of the conventional axial force (x=-A) . 



The component X is positive in the forward direction (thrust) ? 
Y is the side force, positive to the right? Z is the negative 
of the conventional normal force (Z—N^) . The component Z is 

positive in the downwards direction . 

• a?he body axes force coefficients are determined: 



X 

^co S ref 



Y ■ 

q oo S ref 


& 



Z 

g oo S ref 


(60) 


The moment coefficients are likewise defined: 


M 


-<^S£) r c„ = 


My 


X / rex , . c _ i 1 

5 00 S re£ < rsf b ' ^ q^S ref ^ ef 


tef 


■); 


ML 


_ ( jref 


•) 


Cn B q oo S ref' t/ ref ^ 

The lift, drag and side force coefficients are calculated 

relative to stability axes: 

C L = C x sina-C z cosa 


(•61) 


°Y 


= c 


(62) 


(C cosa+C sin a) 

' x z 


and the moment coefficients are obtained: 


C - C cosaH-C sin a ; 

B . 

c = c 

m s mg 

C = - C sincL+c cosa 

n s ^ n B 


(63) 


which completes the linear analysis. 

The above analysis has been programmed, see Appendix B, 
and values for ft for i=l to N appear as output data as well 
as the force and moment coefficients listed above in equations 

(60), (61), (62), and (63). 


•/L 
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Cj. Modified Finite-Step Method: Nonlinear Case 

The analysis of the previous section resulted in a procedure 
for calculating the span loading and integrated wing force and 
moment coefficients which correspond to the linear case. (Actually 
the results were not strictly linear since trigonometric terms 
involving angle of attach and sides lip were not replaced by 'the 
usual small angle approximations) . The linear analysis, however, 
does not predict the stall characteristics of the wing, i.e. it 
leads to an over -estimation of the lift- and an under-estimation 
of the drag at angles of attach above ~10°. It seems reasonable 
to expect that a more realistic estimation of the wing coefficients 
could be achieved by replacing the thin-air foil-theory section char- 
acteristics, (a Q = 2tt l/rad) , with the section aerodynamic character- 
istics as determined from wind tunnel measurements, references 11 or 
12, or as predicted by some of the more advanced airfoil theories 
which tahe into account the viscous forces in the boundary layer, 
see references 13 and 14, While this calculation appears to be 
straightforward, it has not been possible to achieve the desired 
results . 

j The idea of replacing the theoretical airfoil section character 
istics with experimental values is not particularily new, (see 
Reference 12, p. 20) . Sivells and Heely, reference 8, have adopted 
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this procedure using lifting line theory, and it seemed that such 
a procedure would also be adaptable to the Finite-Step Method. 

i 

Ihe procedure adopted for estimating the nonlinear aerody- 


namic coefficients of the wing is basically that of using finite- 
step calculations to estimate the span loading of the wing ini- 
tially. The effective angle of attack at the quarter chord is calcu 
lated from the span loading as in lifting-line theory. These ef- 
fective angles are used to determine the local lift coefficients 
from test data of airfoils. A new span loading can be determined 
using this airfoil lift data and the. effective angles of attack re- 


estimated in an iterative fashion. 

The procedure used in the computer program is amplified below, 
and the difficulties encountered are described. 


The effective angle of attack at each wing segment is 
calculated using the following equation: 


(v di /vj . 

tan a. = • — S3 -Si 

3 •- ' - A 


<W V »> . - 4J 

where the values of are determined from the linear 

3 

analysis. This equation is similar to the effective 
angle of attack calculated in Prandtl’s lifting line 
theory: 


(64) 


a 


'eff . a i 


georn-. 


- w./V 


original RAGE I& 
OF' POOR QUALITY 


(65) 


The airfoil section data c^, c^, c m = f (a) are determined 
at each segment using the value of cu determined from 
equation 64. Data for the tip and root airfoil sections 


are interpolated linearly from root to tip to account 

for the spanwise variation of airfoil section. . The span- 

wise distribution, of circulation is computed from these 

local lift coefficients using: 
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■^original ^ " G -t. °j //2 

4 J 


Values. of r original appear as output from the computer 

r-» 

program. It was intended that the values of T as 
computed in this step he used to replace the linear 
finite-step values used in (1) and these steps repeated 
until no further significant changes in the loading oc- 
curred. However , values of r or ^ g near the wing tips in- 
duced numerical instabilities which necessitated the ad- 
dition steps which are described next. A fuller explana- 
tion of the difficulties is given below. 

3. The values of r^ rig ., as computed' in step 2 above for 

those segments which are located in a region of approximately 
one chord inboard from each of the wing tips, are 
discarded. The finite-step method is then employed to 
calculate replacement values. The procedure used is 
similar to that used in the linear finite-step analysis 
except that the values of the circulation in the central 
portion of the wing are held fixed and equal to ^ or igi na q* 
This is equivalent to replacing lifting-line theory values 
near the wing tips with lifting surface theory values. The 

• ’i 

results of this calculation are labelled L T ,, n and also ap- 

SIMQ * 

pear in the computer output. 
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4 . 


It was necessary to "fair" the values of r SIM Q calculated 
in the tip regions in step 3 -into the central wing values 
determined in step 2 and a polynomial curve fitting 
technique was employed to accompli. sh this. The central 
wing values of ^ SIM q' since they were determined from 
lifting -line theory which is hnown to overestimate the 
lift, were reduced-, using Jones 1 edge correction factor 
obtained from reference 15. This correction is: 


r * 

« _ -orig 

1 SIMQ “ E 


(~ | + c-) < y < (| - c T ) 


where E is the complete elliptical integral with modulus 


3c, i.e.: 

tt/2 

E as V 1 - 1c s sin 3 '§ dq> E(3 c,tt/2) 2 : l 

J o 


and- 


k = 4 1 - 


After application of the Jones 1 factor, the level of T, 
from the finite-step calculations, joined the central 
wing values and the overall distribution could be faired 
using the subroutine MARIAN. The faired circulation 
values were labelled va l ues also appear in 

the computer output. . 
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These stratagems were partially successful in overcoming the 
inherent instabilities which exist in the nonlinear calculation « 

The first instability was uncovered at' low angles of attach 
where the lift coefficient angle of attach relationship is linear. 

The net downwash at the quarter chord of any wing segment, the 
term. w./V in equation (65) , can be regarded as the algebraic 

j m 

sum of the downwash produced by the segment on itself plus the 
upwash contributed by every other segment, assuming they are all 
lifting segments. Suppose that, for some reason, the circulation 
at one particular wing segment is unduly large compared with its 
neighboring segments. When the effective angle of attach is 
calculated on the next iteration, this segment will have a down- 

wash contribution on itself which is too large and « eff will be 
too low. The segments which are adjacent to this segment will exper- 
ience a large upwash and their effective angles of attach and lift co- 
efficients will be too large. As a result, it can be seen that the 
high lift on the first segment will start an oscillation of effective 
angle of attach with attendant changes in section lift and circulation 
upon itself, and that this oscillation will spread laterally to the 
adjacent segments. This oscillation will continue to propagate span- 
wise with each new iteration. This first instability arose because 
of the difficulty of calculating downwash near the wing tips as 
described in step 2 above. 
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Prof. J. Werner suggested that it would be possible to suppress 
this oscillation by using the finite-step method to calculate the 
wing tip loading while retaining the use of section data to deter- 
mine the loading on the central portion of the wing. The polynom- 
ial approximation used to "fair" the loading in the tip and central 
regions served to eliminate the instability. 

The second instability occurred at the higher angles of attack, 
once the airfoil section lift slope became negative. The second 
instability can be visualized by assuming that one wing segment, 
presumably the one with the highest effective angle of attack and 
the lowest stall angle on the wing, reaches an angle where its sec- 
tion lift curve slope is negative. On the next iteration the effec- 
tive angle of attack of this segment will be increased because the 
downwash due its own circulation will decrease, which in turn will 
produce a decrease in lift coefficient. The process is unstable 
and continues until the lift curve slope of this segment becomes 
positive again above the stall angle of attack. The work described 
in references 9 and 10 shows similar results in the steady flow 
calculation. It is indicated in these references that is necessary 
to go to unsteady flow theory in order to achieve even the qualita- 
tive behavior of the wing lift in the stalled region. 

Figure 8 presents the variation of wing lift coefficient with 
angle of attack for an unswept wing of aspect ratio 8.04. The ex- 
perimentally determined pitch characteristic of this wing, compared 
with the results of nonlinear lifting-line theory, are presented in 
reference 16. Both the linear and nonlinear finite— step method 
computer results for five iterations are also shown. The linear 
results as expected, are seen to over-estimate the lift coefficient 
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Since 


by a substantial margin above an angle of attack of 12°. 
the linear cnrve values are used in the first step of the non- 
linear calculation, and they will clearly lead to downwash angles 
Which are too large,- curve #1 is seen to fall below the test data 
by a substantial margin. On the second iteration, the under-estr- 
mated local lift coefficients can result in high effective angles 
cf attack some of which, at the higher geometric angles of attack, 
exceed the local stall angles. Thereafter the stall propagates 
across the wing with the result that the lift coefficients are 
larger than the test data below above a = 14° (they approach the 
linear theory results), but rapidly drop at the higher angles. The 
lift coefficient on the second iteration at a = 22° is off scale 


increases. 



IV. DISCUSSION OF RESULTS 


A parametric study was conducted with the program described 
in the previous sections in order 'to validate the accuracy of 
the methods used. The lift, lift curve slope, and longitudinal 
and directional stability derivatives of a series of wings were 
computed. The aspect ratios considered varied between A = 2 and 
A - 7. Sweepback angles of A = 0 and 45° and taper ratios of 
1 = 0.25, 0.50 and 1.00 were used. Computer runs were made for 
only select combinations of these variables in order to make di- 
rect comparisons with the theoretical and/or experimental results 
available. 

Figures 9, 10 and 11 compare the present results for the 
computation of lift curve slope, static margin, and induced drag 
factor with the vortex-lattice results presented in reference 1. 

It can be seen from these comparisons that the results obtained 
compare quite well with the vortex lattice method results . The re- 
sults should compare exactly for the vortex lattice method using 
one chordwise vortex, (N c = 1) , except for the fact that the vortex 
lattice method uses streamwise vortices in the Treffetz plane and 
moves the "no-flow" boundary condition to this plane also. The 
present method places the vortices on the wing plane. Also, the 
downstream trailing legs lie one-half way between the wing chord 
plane and the free stream velocity direction, as recommended by 
Ruppert, reference 2. The largest discrepancy occurs in Figure 11 
for the swept wings of taper ratio equal to one-half, but percentage 
wise the difference is small. 
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Figures 12 ana 13 present the variation of the pitch damping 
derivatives Oj, and C m with aspect ratio for two sweep angles, 

0 and 45°, tapSr ratios,- X = 1.0 and 0.5*. The computer 

results are compared with the theoretical results obtained from 
reference 17. The agreement between these methods is satisfactory 
except for some differences noted fop the high aspect ratio swept 

wings . 

Figures 14 to 17 present the static later al-directxonal 

. , • , n ssnri n estimated hv the computer 

stability derivatives, , C n and C y ^ , esuimau * 

program with the theoretical and experimental results available 
from various HftCa/HftSa reports, references 7, 18 and 19. 

Figure 14 presents the variation of c. /C L with aspect ratio. 
The results obtained using the present method have been compared 
with three different analytical methods, see Figures 14a and 14b. 

The present method agrees best with the method of reference 7 for 
the unswept wings, which is not surprising since the computer pro- 
gram is based upon the vortex modelling system suggested therein. 
The calculations for the swept and tapered wings, see Figure 14b, 
agree more closely with the results of reference 19, however. 

Figure 14c presents the variation of /C^ with taper ratio for 
the swept and unswept wings of aspect rltio 7, which illustrates 
the better agreement between theory and computer results for the 
highly swept wings. 

A detailed examination of the computer output revealed the ob- 


W standard NASA stability derivatives are used in this section 
of the report and in the figures. 
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vious: that, the rolling moments contributed by the vortices lo- 

cated near the wing tips were the most significant and that the 
integrated values of the rolling moment was the difference in 
the bending moments at the wing root of the left- and right-hand 

wing panels. Thus the value of C. /C T calculated for the straight 

. . . L 

wxng xs equal to the difference of two large numbers of nearly 

identical magnitude. The actual value of /C L for large aspect 

0 

ratio straight wings therefore is sensitive to the minor asymmetry 
of the span loading due to sideslip angle. In the case of the low 
aspect ratio and/or swept wing the differences become more pro- 
nounced and hence better agreement with theory is achieved. 

Separation of the flow from the leeward wing tip due to the 
greater amount of spanwise flow of the boundary layer might cause 
the tip chordwise vortex segment to leave the wing panel and trail 
streamwise rather than chordwise as assumed in the present calcu- 
lations. This effect would produce a significant loading distrib- 
ution change at the leeward tip and alter both the rolling and 
yawing moments due to sideslip angle appreciably. This effect has 
not been investigated numerically. 

Figure 15 presents a comparison of the results obtained with 
the present method with those from references 18 and 20 for esti- 
mating the effect of wing dihedral angle upon C, . It can be seen 

0 

that better agreement is reached with the data and theory (Weiss- 
inger) presented in reference 20. 

Figure 16 presents the estimated directional stability deriva- 
tives for the wings compared with the values obtained from refer- 
ence 18. As pointed out from this reference, the directional 
stability derivative is very difficult to predict and so the dis- 
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agreement shown between the two sources is not unexpected. (It 
should also be pointed out that there is a much more recent and 
advanced source of information concerning the estimation of 
stability derivatives , namely reference 22, and that it is real- 
ized that all of the derivatives presented herein should be com- 
pared with this source. However the methods of reference 22 are 
necessarily longer and more tedious and the line had to be drawn 

somewhere) . 

Figure 17 presents the estimated values of C with those 
obtained using reference 17. The zero sweepbacb values compare 
fairly well but there is considerable disagreement between the 
swept wing values. Reference IS suggests that the theory- of ref- 
erence- 17 is probably not too suitable for estimating of the 

wing. The possibility remains that the present method values are 

reasonable. 

Figure 18 presents the damp-in-roll derivative, C^. Com- 
puted values are shown for a = 0° and 10°. It is seen that the 
values of c. increase slightly with angle of attach. This effect 
■ is not taken into account in the method presented in reference 18 

and deserves further verification. 

Figure 19 presents the variation of the yaw-due-to roll 
derivative, (AC n ) % /C L which contains only the contribution due to 
wing lift and induced drag forces. There is another contribution 
due to the variation of profile drag with angle of attach which of 
course is not taken into account by the present method, linear 
case. Reference 18 indicates that the theoretical values presented 
therein are in agreement with the experimental data, all taper 
ratios, which suggests that further work is required on the present 
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method . 


Figure 20 presents the variation of the side force-due-to 
• roll derivative and fairly good agreement is shown. On the 
basis of the results shown in Figure 19 this agreement might 
he fortuitous, and further study may he warranted here too. 

Figures 21, 22, and 23 present the yaw rate derivatives 
2 2 

/C^, (AC^ ) i a nd C y • Figure 21 shows that the roll-due- 

to~yaw derivative has a variation with aspect ratio which runs count- 
er to the theoretical values given in reference IS. However, as in- 
dicated in this reference, certain empirical corrections are required 
to adjust the theoretical values of C. /C T in order to agree with 

jT ^ 

wind tunnel measured values. On this basis the present method results 
remain questionable and further verification is required. 

The yaw damping derivative, (Ac } /c. 2 , presented in Figure 22 
has the proper variation with aspect ratio but the agreement with the 
level of the theoretical curves is poor for the straight wings and is 
worse for the swept wings. Further investigation here is suggested 
also. 

Figure 23 presents the side force- due- to -yaw derivative and 
good agreement is shown for the straight wing cases while the swept 
wing cases differ by several orders of magnitude. Reference 18 in- 
dicates that experimental data show that the theory (obtained from 
reference 17) is inadequate.. Thus the "good* 1 agreement shown for 
the straight wing cases is suspect. 
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v. CONCLUSION 

The linear finite- step method results for estimating the 
force and moment coefficients in symmetric, sideslipping, and 
rotary flight has been shown to predict values which are, for 
the most part, reasonable. However, it would appear that further 
refinement in the estimation of the lateral-directional rotary 
derivatives is required. 

The. nonlinear finite— step method at low angles of attach 
converges to a solution which is reasonable but does not agree too 
well with the experimental data. At angles of attach above the 
stall the steady flow numerical procedure becomes unstable and is 
unusable in its present form. 
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APPENDIX A: Calculation of the Induced Velocity Influence Coefficents 

The velocity induced at a general point, P, by a horse- 
. shoe vortex having a bound segment with both sweepback and dihedral, 
and trailing legs which are cranked, is developed in this appendix. 
The planview and a sectional view of a horseshoe vortex is shown 
in Figure A~1 * included at the end of this appendix. 

The bound segment, be, lies on the line of quarter chords 
of the appropriate wing panel, see Figure A-l. The bound segment . 
thus has the sweepback and dihedral angles of this line, and is dis- 
placed vertically from the xy reference plane for a wing with non- 
zero dihedral. The chordwise segments, ab and cd, are parallel 
to the x axis. The trailing segments, “^a, and d“ 2 are inclined 
at the angle |3 with respect to the plane of symmetry and their 
projection in the plane of symmetry is inclined at the angle k^a 
with respect to the x axis. For k.^ equal to unity the trailing 
segments are parallel to the free stream velocity vector, while 
for k^ equal to zero these vortices lie in the chord plane of the 
wing with the twist angle removed. The structure of these vortices 
is described further below. 

The wing twist angle has not been simulated by the 
horseshoe vortices, because if it were, the chordwise segments 
of a typical horseshoe vortex would not be parallel for a wing 
with twist. . What has been done is to assume that all the bound 
and chordwise segments on each side of the wing lie in a common 
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plane which contains the quarter chord line and the root chord. 
These planes (one on the left and another on the right sides of 
the wings) are inclined at the angles r with respect to the xy 
reference plane. In order to simulate twist more accurately it 
would have been necessary to introduce a twist angle rotation 
of the chordwise segments and this was considered an unnecessary 
complication. However the wing twist angle is accounted for 
properly when calculating the "no flow" condition at the control 
points . 

th 

The velocity induced at a general point P by the IC 

rectilinear vortex segment, with strength I\ , is derived with the 

aid of Figure A-2. The vectors R, and R are displacement vec- 

K . K th 

tors from the beginning and end, respectively, of the K segment 
to the point p. The unit vector r„ serves to define the dir- 
ection cosines of the segment and is positive according to the 
right hand rule for the assumed direction of The magnitude 

of the velocity induced at P is given by: 

r. 


i 

a K‘ 4hh 


= 


cos 9. + cos 9 

L J 

K \ K K } 


(A-l) . 


where 


h = |r x r | = |r | sin 0 

K K K K 

(A-2) . 

cos 9. = R. * r /Ir, 1 

* l K k k k 

(A— 3) . 
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(A-4) 


cos 0 « ~R_ * r /Ir I 

K K K K 


The velocity at P is calculated: 

«K = % 


.(a-5 ) 


where £ x R r 0 x R 

A _ 2 K K . "K K 

< 3tr — I & . -u- b I — Jl 


|9 • x R 1 
K 


K 


r 2 X ^3 
Z K K 


K 


(A-6) 


hence 


% = In^ cos 9 i k + cos 0 


_ ,,r„ x R 

V V 2 K . 3 k/ 


(A-7) 


The u__, v , w components of q v are found: 
K K K X 


A 


— A 


— A 


“k = V 1 ' ' V K - %‘ 3 ' W K " V k (a - 8 > 

The vectors R, , , and R_, are defined in terms of 

their components 


1 * “2 
K K 


K 


R-i “ R-i t R-t t R-i 
K 1 K 2 k 3 k 


r„ a R • , R„ 


R. 


'2 2 ' *’2 * " 2 
K X K 2 K 3 K 


K = 1,2, 3, 4, 5 


(A-9) 


R_ = R. 
3 K 3 1 


R_ , R_ 
3 3. 

K K 


J K } 


The unit vectors £ can he expresses in terms of the geometry 

K 

presented in Figure A-l. Their components are given in the fol- 


lowing table: 



Where; 


¥ = A, 5=5?., for y.<0 . : : • 

¥ = ~A, I _ _r for y.>0 

♦ 

The vector R for a particular value of K becomes the vector 

J 

r for the next higher value of K, e.g., R-. = ' ant ^ this 

1 j K=2 K=3 

simplifies the calculation somewhat. There are basically only 

six vectors which need to be considered and these are the vectors 

which mark the ends or corners of the cranked horseshoe vortices, 

namely: R , R^ R^, R,, Rj, and 1^. It is noted that the ends 

of the trailing legs are not taken to be located at infinity but 

rather at approximately 1000 root chords downstream for ease of 

computation. The value cos 0^ for K^, and cos 0^ for K=5, are 

nearly equal to unity at such a great distance downstream. The 

aerodynamic influence of the trailing legs at a distance beyond 

1000 c is negligible. The components of R and R for K=l— -5 
R K 

are given in Table A-2. Notice that x,y,z have been replaced with 
•f,g,h, respectively in this table since the latter variables 
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will be the ones used in the formulation given in the main body 
of the report, see equation 24. 


Table A— 2 Components of Rx K and Ra K for K -l-*5 


K 

» 1 
H 

N 

1 

i 

,R3 k 

Rx 

Rx 

2 k 

Rx 

__±K_ 

Ra 

1 K 


Rs a 

K 

K=i 

1000c cosp* 

XV 

cosOcx a ) +s* 

cos<|)tan¥ -f 
y * C MI +f 

1000c sing 

XV 

4- s coslj) 

+g 

1000c cosp 
i‘sin (kj a) 4 
s sin^) 

+ h 

s cos(i)tan¥ 

. MI 
+ f 

s cps$ + g 

s sin(j) + h 

K=2 

1 

= \ 
X 1 

\ 

1 

R 

3 1 

s cos (j) tan Y 
+ f 

s cos^) 4- g 

s sin(|) + h 

K=3 

R-3 

■ . ^ ' 

* 

\ 

2 

\ . 

-s’ cos^tanY 
+ .£ 

s cosljl + g 

-s sin(j) 
4- h 

K=4 

R 3 

l 3 

R 3 

2 3 

R 3 

, 3 3 

-s cosi* 
tan! + 
^ C PL 4- f 

— s cos(j) 

+ g 

-s sin(|) +h 

.V 

K=5 ' 

R 3 

V 

; -i 

R3 2 4 


1000 c cosp 

XV 

-cosfc^cx) - 

s cos^tanY 
* ?*C PL + * 

1000c sin |3 

■ R 

— s cos® 4-g 

' 

1000c cosB 
R 

-sin (k^ a) - 
s sin(j) + h 


■i. ; , . 

tv _ The induced velocity components given by equation A-8 

i 

are found by performing the steps indicated in equations A-2, A-3, 

A-4, and A-7„ For the sake of completeness the amplified equa- 
• • 

» 

tions used are given below. 
























"■«- *», - ■>, v'* ' w- w* ww 


“O ~'Q 3 2 

K K K K 


R 1 ^ R 1 ^ ~ 

K K 


R «'|» 1 = 

K 3 K 


R 3 + K 3 -F R* 
x i 3 • 

K K K 


R 3 + R 3 + R- 3 

J 2 ^3 

X K K K 


K K K K 
~)h 
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K K K K 
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cos 0 1 « 
X K 


R, R , + R, R 0 + B B 

X K X K K K K K 


/B-, 


K 


(A— 13) 


cos 0_ = ~ 

3 k 


R n R, + B R + R R 

1 2 1 2 2 3 3 

K K K K K K 


/B, 


K 


(A-14) 


U K 4TT’ 


r. /cos 01 + COS 0a J\ 

ix ■ 
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(R R - R R ) 
Z 2 3 2 

j k k k 


(A-15) 


r . ft COS 01 „ + COS 0a . 


V K 4tt\ 


K 


(E 2 R 3 ‘ E 2 E 3 } 
J K K K K 


(A~16) 


r. cos 0i + cos 0a « 

W K= iW V (R 2 V - R R 3 ) 

K K K K K 


(A-17) 



(A-18) 


(A-19) 


(A-20); 


! 

I 

{ A~10 
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The influence coefficents are defined: 


r. 

u. = tt: Fu. . 
1 4TT XJ 

r. 

V. - T” Fv. . 
3 4TT X3 


(A-23 ) 


(A-22 ) 


W. = -rrr Fw. . 

3 4TT 13 


(A-23) 


hence 



Fv. .= 
3-3 


Fw. .= 
2-3 



K =1 

5 



cos 

01 K + 

cos 

a 

CO 

CD 


K 



cos 

01 K + 

cos 

03 K 


h K 



cos 

CD 

+ 

cos 

03 K 


, -Tv- - ■C'-o 

' 2 k \ \ 2JU-24) 


K 


'3 1 

J K K 


'1 2 
K K 


\ 3 3j(A-25) 


>9 

2 k V 1 *- 261 


K=1 


In the calculation of the influence coefficents, equations A-24, 

A-25, and A-26, for the case where the general point P corresponds 
to the point d^ (which lies on the bound vortex segment, see Sec- 
tion III.B.3.b.) it is necessary to omit the terms for K=3 when 
i=j a In this case and others where the numerical evaluation 
of the influence coefficents may have become difficult owing to 
^ 2^0 fg_ct that 0), -» 0, x.e. , the point P approaches .the lxne upon which 
the segment lies and equation A-l becomes indeterminate, the 
computer omits the induced velocity of that segment. A general 
test is made and the computer omits the induced velocity com- 
ponents of any vortex segment when the value h K becomes less 
than 0 . 1 s. 
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FIG, A— I PLAN AND SECTIONAL VIEW OF A 

TYPICAL HORSESHOE VORTEX 
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1. Program Description 

The program' computes the spanwise distribution of circulation. 
r(y) “ r./V , and the linear force and moment coefficients for 
a wing alone in a steady, incompressible inviscid fluid flow using 
the finite-step method- Provisions are included for studying 
variations of angle of attack, angle of sideslip, and the angular 

rates of rotation P, Q, and R. 

The wing geometry is limited to a single trapezoidal panel 

per side and left-to-right geometric symmetry is assumed, except 
that differential wing twist is considered in order to permit the 
study of manufacturing anomalies. The wing is divided into an 
even number, N, of segments of equal spanwise dimensions. The 
number chosen for H can be varied between twenty and one hundred 
in the program. The segments are numbered 1 to N, from left to 
right, beginning on the left wing tip. The program solves N 
equations for the spanwise distribution of circulation in every 
case, advantage could have been taken of aerodynamic symmetry 
for those cases dealing with symmetric flight only, and H/2 
equations solved instead. The program is inefficient to this 

extent. 

A provision has also been made for the calculation of the 
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nonlinear aerodynamic coefficients of the wing but this part of 
the program is not working satisfactorily as described previously. 
This calculation may be deleted as described in section 3 below. 
There are also some minor variations provided for preparing the 
input data concerned with wing geometry . 

2. Operating Information 

Core and Time Requirements: 

Computer: CDC 6600 

Core : 77 K s to load 

70 K_ to execute 

s 

Time : Approximately 4 minutes 

3. Input Data 

There are two options for loading the basic wing geometry. 

In the first option, the span and root and tip chord dimensions 
are given and the planform area, aspect ratio, and taper ratio 
are computed. In the second option, the area, aspect ratio and 
taper ratio are loaded directly. The root chord is assumed equal 

to unity in the second option. 

The angles of attack and sideslip and the angular rates P, 

Q, R are varied in the following manner. Five D0 L00PS are nested 
from the outermost variable R, to the innermost variable j3, xn the 
order R, Q, P , a, 0. The calculations proceed by selecting the 


62 



first set of values, namely, Q,- V *" a varies 6 fr ° m B 1 fc ° 

B ib . Upon completion of this calculation, the next a is chosen 
holding Q r fj. «*- and calculations are made for the various 

0 values. After calculations have been made for all the angle 
of attack values considered, the computer then repeats these 
calculations for the different roll rates. Subsequently, all 
the preceding calculations are repeated for the different pitch 
■ rates and then the different yaw rates in that order. It is more 
efficient usually to use different cases to consider separately 
the effects of these variables than to use the nesting provided 

in the program. 

Multiple cases can be run consecutively by repeating the 
information retired for a single case, see description given for 

card ® below „ 


Cell? cl 

No, Variable Format 

• © N0CSSE 15 


© CRF 


8F10 . 0 


Description . 

Number of cases to be run. The 
program is presently arranged so that 

only cards @ to © or @ can be 
modified. Changes to cards © , etc., 

requires a complete reloading of the in- 

/ 

put. 

c „ mean aerodynamic chord, or 
ref 
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Card 

No. 


Variable Format 



( Contd . ) 

CLD 

CGD 

ETLD 

ETKD 

D340D 


8F10 . 0 


tl 


II 


It 



Dihedral angle of quarter chord 

line, T, degrees, positive for tip 

chord elevated above root chord. 

Twist angle of the left wing tip 

chord, e , degrees, positive for 
R 

washin (negative for washout;) . 

Twist angle of the right wing tip 
chord, e_ , degrees, positive for 

i R 

washin. 

The angle the mean camber line mahes 
with the chord line at the wing root 
(y-0) , ( 6 3 / 4 ) 0 ' degrees, positive 

for positive camber, see Figure 4. 
This value is equal to the geometric 
value for the linear case only, and 
should equal the negative of (o^ ) Q ' 
approximately when the airfoil aero- 
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Variable Format 


Description 


D34TD 


XCG 


YCG 


ZCG 


SKI 


8F10 . 0 


tt 


11 


dynamic data is known. 

The same as D340D , except that this 
value pertains to the wing tip rather 
than the wing root. 

The longitudinal position of the 
center of gravity, measured from. the 
quarter chord point of the root chord, 
x CG , (m/ft.) , positive for the CG 
ahead of the \ C point, see Figure 1. 
The lateral displacement of the CG 
from the plane of symmetry, y , 

(m/ft.) , positive for a displacement 
to the right. 

The vertical displacement of the CG 
from the x,y reference plane, z , 
(m/ft.), positive for a displacement 
below the reference plane, see Figure 2. 
The value , k , where k is the 

X j. 

fractional value of the displacement 
angle of the wing wake above the wing 
chord plane; when k^O, the wake trails 


ORIGINAL, Maw k . 



Card 

No. 


Variable 


Format 


© ( Contd . ) 


t 

4" 


ALOKD 


8F10 . 0 


ALOTD 


N 


ITA 


1015 


Las 


* ' 

' oS * 


•Description 


in the -wing chord plane and when 
the wake trails parallel to 
V „ Nominal value is k -~h- 


CO 


angle of zero lift of the root chord, 

(a ) r degrees, obtained from airfoil 

b 0 0 

section data. Value should be negative 

for a positive camber. 

Angle of zero lift of the tip chord, 

(cl ) r degrees, from airfoil section 
"0 0 

data. 

Number of wing segments, 20 < N < 100, 
^nominal" 

Number of angles of attack used m 
the table of airfoil section character- 
istics, see card ® below. The 
maximum value for IT A is 20. 

Option for loading planform geometry, 
LAS <* 0 calls for card ® whereon 
span, root and tip chord, and arr 
speed are placed. For LAS > 0, card 
( 5 ^ is used instead, wheron reference 
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Card 

No. 


Variable Format 

© (contd. ) 


IA 1015 


IB 


IP 



> 


'■IQ 

It 

IR 

> 

LIN 

It 



Description 

area, aspect ratio , taper ratio and 

air speed are placed. 

The number of angles of attach to 

be considered are indicated by the 

variable IA, IA “25. 

max 

The number of angles of sideslip 

to be considered are indicated by the 

variable IB, IB =5. 

max 

The number of roll rates, pitch 
rates, and yaw rates to be considered 
are indicated by the variables IP, 

IQ and IR, respectively. A maximum 
of five values of each may. be considered. 
The variable LIN exercises the option 
to calculate the linear aerodynamic 
results only, LIN = +1.7 or the combined 
linear and nonlinear results, LIN = -1. 

In the linear only case, it is un- 
necessary to lead cards (s) to (25) . 


The following card is an optional card, but either it or 
card (5^) must be used. In order for (5a) to be used. 
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< 0, see card (4) . 


Variable Format 


Description 


8F10.0 Wing span, b, (m/ft.) 


Tip chord, C T , (m/ft.) 

Root chord, C R , (m/ft.) 

Air speed, V , (m/sec. or ft. /sec.) 
Wing planform area, , (m 3 /ft 3 ) 

Aspect ratio/ A, nondim. 

Taper ratio, X, nondim. 

Air speed, V , (m/sec. or ft. /sec.) 


The airfoil section characteristics are placed on the next 
group of cards which are ITA in number. The root chord 
values are listed first, followed by the tip values. If 
LIN >0, omit these cards. 


TA(I) 

CLR(I) 


CDR(I) 


7F10.4 Angle of attach, a, degrees 


Root section lift coefficient, c^, 
nondim. 

Root section drag coefficient, c^, 
nondim. 
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Card 

No. 




Variable 


Description 

CMR (I) 

7F10.4 

Root section pitching moment co- 

CLT(I) 

II 

efficient, c . , nondim. 

ra c /4 

Tip section lift coefficient, c^ , 



nondim. ' ’ 

CDT (I) 

It 

Tip section drag coefficient, c^ , 



nondim. 

CMT(I) 

II 

Tip section pitching moment coefficient. 



c • , nondim. 

ra c/4 


The next group of cards list the values of the angle % s ) of 
attach, angle (s) of sideslip, roll rates, pitch rates, and 
yaw rates for which calculations are to he made. 


@ aii(i) 


if 

@ AL (IA) 


8F10.0 Value (s) of the angle (s) of attack, a, 
degrees for which calculations are 
to he made. A maximum of 25 values 
may he considered. 
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Card 

No. 



© 

0 


Variable 
BE (1) 

i 


A 


be Tib) 


PCX) 

T 

P(IP) 

■*% 


Q(l) 

Y % 

Q(IQ) 

R(l) 

R(IR) 


Format 

8F10.0 


It 


tl 


II 


Description . — * 

Value (s) of the angle (s) of sideslip, 
p, degrees. A. maximum of 5 values 
may he considered. 

Value (s) of the roll rate(s) , P, rads. /sec 
A maximum of 5 values may he considered. 
Value (s) of the pitch rate(s), Q, 
rads. /sec. 

A maximum of 5 values may he considered. 
Value (s) of the yaw rate (s) , R, 
rad. /sec. 

A maximum of 5 values may he considered. 


Additional cases may he 


listed hy repeating the input cards listed 


above from 2 to 5a or 5h . 
arrangement of the input cards. 


Table , B-l presents the general 


4. Output Data 

The output will he described for the case where the combined 
linear/nonlinear option has been chosen. (In linear case only, 
the nonlinear output will, of course, not appear.) 

The output is arranged in four basic groupings: 

a. The wing geometric variables are listed together with 


some input constants which serve to identify the case: 

CREF « c ref = c, ETL = e° , (DELTA 3/4)0 = (6^^) g , 

XCG a x cG/ ZCG = z cq/ GAMMA = T ° , ETR a e° R , (DELTA 3/4) T 

(5 3/4>T ' YCG = Y CG' n = V = A ° = A °/4 ' 

ALPHA = a° / V INF = V (ft. /sec. or m/sec.), CT = c , 

■LiO 03 T 

B = b, CR = c (or alternatively SREF = S _, A = A, , 

K. ref 

LAMBDA = X) , E = E. These values are not repeated any- 
where in the run. 

The angular velocity components R, Q, P (rads. /sec.) 
and the angle of attack and angle of sideslip are listed 
on the next line. 

The linear data is listed in the following manner: 

1- The linear-analysis spanwise values of the circulation 
r\ = (T /VJ is listed left to right, for the N 
segments, beginning of course with the value 1=1, 

continuing to. I = N. (1 = 1 corresponds to the segment 
on left wing tip and I — N corresponds to the segment 
on the right wing tip.) 

2. The linear, integrated, force coefficients c , C 

x y 

C and the moment coefficients: CLB = C , CMB = c 

■Cb m B ‘ 

and CHB = C for the body axes are printed on the 
next line. A title appears in the preceding line. 

3. The above data is transferred to stability axes. 
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a. 


labelled as such, and printed on a single line. 

The force coefficients are CL = C , CY = C . , arid 

L y 

CD = C ■ the moment coefficients are CLS = C , 

induced 

CMS = C , and CMS = c . 
m s n s 

The nonlinear data is listed in the following manner: 

/vf 1 

1. The M values of r original are printed for the wing 
segments, beginning with the left tip value and 
ending with the right tip value, from left to right 
across the page. These values correspond to the 
first iteration. 

2« In a similar fashion, the values T . are printed. 

Sxmq ' 

followed by the values of f 

Mar ran 

3. The N values of the pairs of number TCLJ(J) == 2T./c. 

^ 3 

and CLJ(J) = c are printed side-by-side. The vhlue 

TCLJ(J) is the "theoretical" section lift coefficient 
. tii 

at the j segment corresponding to the value of the 
circulation (at thht' segment) determined by linear 
theory, while CLJ(J) is the section lift coefficient 
at the j segment determined from the wind tunnel 
data tabulations corresponding to the computed value 

of « . . 

3 

4. The effective angle of attack, aj , is printed for 


ORIGINAL BAGS 1$ 
OF POOR QUALITY 
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each of the N segments. 

5. The integrated nonlinear force and momeht coefficients 
for the body and stability axes, as was done for' the 


linear calculation, are printed next. 

6. Four more iterations follow the first before the 
angle of attach or angle of sideslip change {or 
P, Q, R change) to the next value. 

5 . Programming information 

Purpose of Subroutines: 

EEL2 Computes the complete elliptical integral 

E(h, tt/ 2) where k=^/l-{b/a) 2 and a/b=Tr A/4, for 
correcting lifting- line theory according to 


R.T. Jones: c ^ =2TTa e ff //E * 

SIMQ Solution of simultaneous linear algebraic equations 

used primarily for the determination of the r vector 
which satisfies the no-flow boundary condition in the 
linear calculation, but also used in the nonlinear 


calculation for determining the lifting surface value 
of the circulation in a region one tip chord in length 
located at the wing tips. 

FUVW Calculates the induced-flow velocity influence coef- 

ficients . 

MARIAN A seventh order polynomial curve-fitting procedure 

which smooths the spanwise distribution of circulation. 






0 






& 


There 


joining the section values in -board to the lifting- 
surface theory values near the tips, 
is no interdependence of these subroutines. 


73 



»1 ' 
II . 

A & 

& 'S' 

. IV 


•- • TABLE B l- . ,, 


(7) "A VOtSS^ *• » H »-«UQJU-33. 

^ ^ L_ ! — J 1 — i — * *— 1 — > I 1 t i l I i t i . . . . I 


JU2JJUUUUJ.ilJJJ.<J 


»-». 7 * 7,' JJ.1 


» T4 71 TI7TM7 



75 



FIG. 2 LEFT HAND WING PANEL DEFINING 
TWIST ANGLE {ZERO DIHEDRAL 
CASE) 
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FIG. 3 LEFT HAND WING PANEL DEFINING 
UNIT TANGENT AND NORMAL VECTORS 
ALONG THREE-QUARTER CHORD LINE 
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MEAN CAMBER 
LINE 



z 


FIG. 4 TYPICAL CHORD PLANE DEFINING 
BOUNDARY CONDITION AT THREE - 
QUARTER CHORD POINT 
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FIG. 5 RIGHT HAND WING PANEL DEFINING 

UNIT TANGENT AND NORMAL VECTORS 
ALONG THREE-QUARTER CHORD LINE 


ORIGINAL SAGE IS 
OP POOR QUAIOT 
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FIG. 7 SCHEMATIC DRAWING ILLUSTRATING 
BOUND AND CHORDWISE VORTEX 
SEGMENTS FOR CALCULATING FORCES 
AND MOMENTS 




a , ANGLE OF ATTACK 

FIG. 8 VARIATION OF LIFT COEFFICIENT WITH 
ANGLE OF ATTACK - NONLINEAR THEORY 
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FIG. 9 COMPARISON OF PRESENT 

VORTEX LATTICE METHOD - STRAIGHT WINGS 




ORIGINAL BAGS 15 
OF POOR QUALITY 



FIG 10 COMPARISON OF PRESEMT CALCULATIONS WITH 
VORTEX LATTICE METHOD, SWEPT (A=45°), 

UNTAPERED WINGS 




FIG. 10 CONCLUDED 


B6 



FIG. II COMPARISON OF PRESENT CALCULATIONS WITH 
VORTEX LATTICE METHOD, SWEPT (A = 45°), 
TAPERED ( X »0.5) WINGS 








A, ASPECT RATIO 


FIG. 12 COMPARISON OF PITCH 
C Lq VS A 


2 


A, ASPECT RATIO 


IMPING DERIVATIVES, 





FIG 13 COMPARISON OF PITCH DAMPING DERIVATIVES, 
Cm q VS A 



0 



X =i.oo 
A =o° 



O PRESENT METHOD 


o 1 1 1 

2 4 6 8 


A, ASPECT RATIO 


X *1.00 

A =45° 



A, ASPECT RATIO 


FIG. 14 a COMPARISON OF PRESENT METHOD WITH RESULTS OF 
VARIOUS REFERENCES FOR ESTIMATING C/ /Ci • 
UNTAPERED WINGS P 



A, ASPECT RATIO 


A, ASPECT RATIO 


FIG. 14 b CONTINUED 1 TAPERED WINGS 
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I? 


- 0.2 


.25,, 5 ,1.0 NACA TR 1098 

1.0 • WEISSINGER TR 969 
1.0 ) 

0.5 PRESENT METHOD 
0.25 ’ 

1 . 0 EXPERIMENTAL 
DATA TR 969 


A, ASPECT RATIO 


A, ASPECT RATIO 


FIG. 15 COMPARISON OF VARIOUS METHODS FOR ESTIMATING THE 
EFFECT OF WING DIHEDRAL ON C/^ 




Cn/C L 2 


FIG. 16 DIRECTIONAL STABILITY DERIVATIVE 


0.4' 


A=0° 


Jl 

03- 1.0 NACA TR 1581 

C Y o ° '‘^PRESENT METHOD 



FIG. IT SIDE FORCE - DUE- TO 


8 


2 


4 


6 


0 


A, ASPECT RATIO 


SIDE SLIP DERIVATIVE 



FIG. 18 COMPARISON OF PRESENT METHOD WITH METHOD OF 
NACA TR 1098 FOR ESTIMATING DAMPING - IN- ROLL C/ 





A=o° 


-0.5H 



A, ASPECT RATIO ' A > AS P ECT RATIO 


FIG. 19 YAW- DUE- TO ROLL DERIVATIVE, (AC np )|/C L DUE TO 
LIFT AND INDUCED -DRAG FORCES ONLY, c.q. AT a.c. 





'& 00<1 


X. 

L0 y NACA TR 1098 
0.5 J 

□ as} present method 


O 


□ 



i 

I 




0 


it? 


J 1 L 

2 4 6 

A , ASPECT RATIO 


J 

8 


0 


2 4 6 

A, ASPECT RATIO 


8 


FIG. 21 ROLL— DUE —TO YAW DERIVATIVE 

C L 
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A, ASPECT RATIO A, ASPECT RATIO 


FIG. 22 DAMPING IN YAW DERIVATIVE, (&Cn r ) DUE TO LIFT 
AND INDUCED DRAG FORCES ONLY, c.g. AT a.c. 
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